How to use retarded Green's functions in de Sitter spacetime 



Atsushi HIGUCHI 1 and LEE Yen Cheong 2 

Department of Mathematics, University of York, 
Heslington, York YO10 5DD, United Kingdom 
email: 1 ah28@york.ac.uk, 2 yl538@york.ac.uk 
(Dated: September 18, 2008) 

Abstract 

We demonstrate in examples that the covariant retarded Green's functions in electromagnetism 
and linearized gravity work as expected in de Sitter spacetime. We first clarify how retarded Green's 
functions should be used in spacetimes with spacelike past infinity such as de Sitter spacetime. In 
particular, we remind the reader of a general formula which gives the field for given initial data on a 
Cauchy surface and a given source (a charge or stress-energy tensor distribution) in its future. We 
then apply this formula to three examples: (i) electromagnetism in the future of a Cauchy surface 
in Minkowski spacetime, (ii) electromagnetism in de Sitter spacetime, and (iii) linearized gravity in 
de Sitter spacetime. In each example the field is reproduced correctly as predicted by the general 
argument. In the third example we construct a linearized gravitational field from two equal point 
masses located at the "North and South Poles" which is non-singular on the cosmological horizon 
and satisfies a covariant gauge condition and show that this field is reproduced by the retarded 
Green's function with corresponding gauge parameters. 

PACS numbers: 04.62. +v 
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I. INTRODUCTION 



Over the past few decades quantum field theory in de Sitter spacetime has been devel- 
oped mainly because of its relevance to the inflationary universe scenario [l|, H, H, 0, 0, @] 
(see also Ref. Q). In particular, infrared properties of linearized gravitational field have 
attracted much attention. The mode functions for linearized gravity are similar to those 
of the minimally-coupled massless scalar field if one imposes gauge conditions natural to 
the spatially-flat coordinate system of de Sitter spacetime The two-point function of 
the minimally-coupled massless scalar field is known to be infrared (IR) divergent (see, e.g. 
Refs. BSEjini), and as a result one finds that the two-point function for linearized grav- 
ity is also IR divergent in this gauge. However, since linearized gravity has gauge invariance, 
one needs to determine whether or not these IR divergences are a gauge artifact. Indeed 
it has been shown that they can be gauged away (by gauge transformations which are not 



localized) [12], [13| . It has also been shown that the two-point functions in gauges natural to 



other coordinate systems are IR finite [14|, [15 



Infrared finite graviton Feynman propagators in the covariant gauge have also been con- 
structed [l6|, 17]. However, it has been claimed that the covariant Feynman propagators 
lead to incorrect retarded Green's functions which violate the linearized Einstein equa- 
tions 18, This assertion is partly based on the work of Antoniadis and Mottola [20], 
who claimed that in order for the scalar metric perturbation to obey the equations of motion, 
the de Sitter invariant construction of the retarded Green's function must be abandoned (see 
also Ref. 2lj]). However, it had been pointed out by Allen j22| that the value of one of the 
gauge parameters used in their work introduces spurious IR divergences for the Feynman 
propagator in the sense that the IR divergences are absent for other values of gauge pa- 



rameters, e.g. those adopted in Refs. |16|, ll7|. As far as non-interacting linearized gravity 
is concerned, the conclusion of Antoniadis and Mottola is entirely due to these spurious 
divergences. Moreover, by a closer inspection one finds that their retarded Green's function 
for the scalar metric perturbation, which is IR finite unlike the Feynman propagator, does 
satisfy the correct field equation contrary to their original claim [23j . Thus, any objec- 
tion to the covariant retarded Green's function for linearized gravity based on their work is 
completely unfounded. 

Another feature of the covariant retarded Green's functions that might appear to cast 
doubts on their validity is that, since they are causal, the field generated by an inertial 
point source in de Sitter spacetime has support only in half of the spacetime (l8j . Thus, for 
example, the electric field generated by a freely falling point charge cannot satisfy the Gauss 
law because the flux out of a sufficiently large sphere around the charge vanishes, the field 
itself being zero there. This apparent paradox can be resolved by recalling that in general 
the retarded Green's function is used to generate the field in the future of any Cauchy surface 
for given initial data on it and a source in its future. Thus, in order to reproduce the field 
using the retarded Green's function one needs to include the contribution from the initial 
data on spacelike past infinity in de Sitter spacetime as well as that from the source: the field 
in de Sitter spacetime is not determined by the source alone but is influenced also by the 
initial data on past infinity. There may be philosophical uneasiness about the fact that the 
electric field on past infinity is necessarily nonzero if charged particles are present there j24[ , 
and some authors have proposed that this field should be generated from the source by 
using half-advanced Green's function [25T ]. (A related fact jU 25[ is that an electric charge 
does not generate a purely retarded electromagnetic field with support in its causal future.) 
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However, the retarded Green's functions, as a mathematical tool, are obviously not designed 
to generate the initial data on past infinity. 

Thus, we find no compelling reason to doubt the validity of covariant retarded Green's 
functions either in electromagnetism or linearized gravity in de Sitter spacetime. However, 
in view of the recent claim that they do not work properly, it will be useful to demonstrate 
how they work in some examples. The purpose of this paper is to show that the covariant 
retarded Green's functions for electromagnetism [26] and linearized gravity 17j do indeed 
generate the fields which obey the equations of motion if we use the formula involving not 
only the source but also the initial data on a Cauchy surface. Our calculations also serve as 
a check for the Feynman propagator obtained in Ref. [13]. (See Ref. 27] for a calculation 
of a physical quantity, the Weyl-tensor correlation function, using the two-point function of 
Ref. Q.) 

The rest of the paper is organized as follows. In Sec. Owe address by using Carter- Penrose 
diagrams the necessity of including the contribution from past infinity as well as that from 
the source. Then, we review the derivation of a general formula for reproducing the field 
from the initial data on a Cauchy surface and a source in its future in terms of the retarded 
Green's function. We then apply this formula to three examples. In order to illustrate its use 
in a simple setting, we first discuss how the electric field with a static charge is reproduced 
in the future of the t = Cauchy surface of Minkowski spacetime using the Cauchy data 
on the t = hypersurface in Sec. IIHI In Sec. [IV] we reproduce the electromagnetic field 
with two charges, one at the "North Pole" and the other with the opposite sign at the 
"South Pole" . We demonstrate that the general formula in Sec. [IT] correctly reproduces the 
field throughout the whole manifold of de Sitter spacetime. In Sec. [V] we first write down 
the linearized gravitational field with two equal mass points at the North and South Poles 
smooth at the horizon and satisfying a covariant gauge condition. Then we show that the 
linearized gravitational field on the entire spacetime is recovered using the general formula 
in Sec. II as expected. Throughout this paper we call the contribution to the field coming 
from the source the source field and the contribution coming from the initial data the initial 
field and use natural Planck units h — c—G—1 and the sign convention — h ++. 



II. GENERAL FORMULA FOR THE USE OF RETARDED GREEN'S FUNC- 
TION 

The electromagnetic and gravitational fields generated by the retarded Green's function 
propagate at most at the speed of light. On the other hand the Gauss law in electromag- 
netism and a similar law in linearized gravity in spacetime with a Killing vector field imply 
that the total outward flux of the field out of an arbitrary closed spatial surface is equal 
to the total conserved charge enclosed by such a surface. These two facts would appear 
to invalidate the use of covariant retarded Green's functions in de Sitter spacetime [181 ] . 
This apparent paradox can be best understood using Carter- Penrose diagrams. The Carter- 
Penrose diagrams for Minkowski and de Sitter spacetimes are shown in Fig. [1] (see, e.g. 
Ref. [28l]). The bold line represents the world line of an inertial source originating at past 
infinity i~ (X~) and ending at future infinity i + (X + ) in Minkowski (de Sitter) spacetime. 
The field causally generated by the source along its world line has support in the shaded 
region. For electromagnetism, for example, the field generated using the retarded Green's 
function by the charge cannot satisfy the Gauss law on any constant-time hypersurface in 
de Sitter spacetime since the electromagnetic field vanishes outside a sphere of some radius 
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(i) (ii) 

FIG. 1: (i) The Carter-Penrose diagram of Minkowski spacetime; (ii) Carter-Penrose diagram of 
de Sitter spacetime [see Eq. ()4. 1[) for the definition of x]- The field coming from a classical source 
(its world line represented by a vertical bold line) is nonzero in the shaded region. 



at each time, as can easily be seen by drawing a horizontal line across the Carter-Penrose 
diagram, whereas in Minkowski spacetime the Gauss law is satisfied by the field generated 
by the retarded Green's function from the source. The cause of this apparent difficulty in 
de Sitter spacetime is the spacelike nature of past infinity I~ . 

In fact there is nothing wrong with the retarded Green's functions de Sitter spacetime. 
In some sense this apparent difficulty is caused because one is implicitly requiring too much 
of the retarded Green's function. To understand this point one needs to recall how the 
retarded Green's function is used to generate the field in the future of an arbitrary Cauchy 
surface. There is a formula, which is the key formula of this paper, for generating the field 
in terms of the retarded Green's function for a given initial data on a Cauchy surface and a 
given source. This formula is well known for scalar fields (see, e.g. Ref. 0]) and is known 
for other bosonic fields as well, but it will be useful to review its derivation here. 

Suppose the Lagrangian density for a local field Aj (with an upper case Latin index 
representing a set of indices) is given by 

C = ^ [T aIbJ V a AjV b Aj + S IJ A T Aj] , (2.1) 

where T aIbJ and S IJ are tensors satisfying T aIbJ = T hJaI and S IJ = S JI and independent 
of Aj. The conjugate momentum current is defined by 

1 r)C 

Let L n be the differential operator defined by the equation 

(L^AY 1 = tt c/ . (2.3) 

Thus, the operator L v maps the field Aj to its conjugate momentum current ir cI . The 
Euler-Lagrange equations read 

L IJ Aj = V c vr c/ - S IJ Aj = . (2.4) 
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We assume that there are unique retarded and advanced Green functions, Gf r (x, x') and 
Gfj,(x,x'), satisfying 

Li J Gf/{x, x') = 5 V 4 (x, x') , (2.5) 

where the (un)primed indices refer to the (un)primed point (x) x' and where the differential 
operator iJJ acts at point x. The retarded (advanced) Green's function is also required to 
vanish if point x is not in the causal future (past) of point x', i.e. if there is no future-directed 
(past-directed) causal curve from x' to x. The delta function 5 1 p5 A (x, x') is defined by the 
property that 

J d*x' y /=g i 8 I I ,8 A (x, x')A r (x') = A\x) (2.6) 
if A 1 \x) is smooth and compactly supported. It can be shown, as is well known, that 



I& J 'G?j,{x,a?) = 5 v iP{rf,x) 



I' jt4/„ / 



(2.7) 



which implies Gfj,(x,x') = Gf,j(x',x) by the assumed uniqueness of the advanced Green's 
function. We present a proof of this formula in Appendix [A] 

Let us denote by -D + (S) the future domain of dependence of the Cauchy surface E. (See, 
e.g. Refs. 28, 3(| for the definition of the domain of dependence.) Let x be a point in the 
future of E, i.e. x G _D + (E), and x' be a point on the Cauchy surface E, i.e. x' G E. Consider 
a source represented by the current J 1 that is nonzero in the region D + (E) . The unique 
solution to the inhomogeneous equations of motion 



L IK A K = V c tt c/ - S IK A K = J 1 



(2.8) 



in the future of a Cauchy surface E with given initial data on E can be expressed in terms of 
the retarded Green's function Gfj,(x,x'). Let (Aj/,tt cI N c >) be the initial data on E, where 
N c r is the past-directed unit normal on E. (Note that the time component of N c i is positive.) 
Let E' be another Cauchy surface in the future of E and let x G U = -D + (E) r\D~(T,'), where 
D~(E') is the past domain of dependence of E'. Then 



AAx) 



d i x'y/^(yj^(x,x , )5 J I 'Ap(x') 

= [ d'x'^i^L^'GfAx^^Mx'), (2.9) 
Ju 

where we have used Eq. (|2.7p . Using Eq. (12.81) . we find that Eq. (I2.9p can be expressed as 



Aj(x) = I dVv C 0OcO{V c ' \{L v G R )i {x 1 x')Ar{x , )-Gf JI (xiX , )n c!j '{x') 
+Gf II (x,x')J I '(x') 



(L^ i? ) / c ' / '(x,x')V c ^ / ,(x') - V c , {Gfj,(x,x')) n c ' J '(x') 



(2.10) 



By Eq. (T22D and the equation (L v G R )i {x,x') = T c 1 a J V a 'Gfj,(x, x') together with the 
symmetry property of T c 1 a J , the last two terms on the right-hand side of Eq. (12.101) cancel 
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out. Then, by the generalized Gauss theorem we find 



AAx) 



V 

+ I dY> c < 
i 

dTi r i 



{L^G H )r (x,x')A r (x') -G?j,(x,x')t[ c ' j '(x') 



{L^G R )r {x,x')A P {x') -G?j,{x,x')ir c ' J '{x') 



c'l 



(2.11) 



where d£ c / = dY,N c >. The second term, the integral over vanishes because Gfj,(x, x') = 
if x is not in the causal future of x' . Hence, 



Aj(x) = A i I s \x) + 4 I \x), 
where the source field A^p (x) and the initial field A^P (x) are given by 



Af\x) 
A { /\x) 



d i x'^/-g{x')Gf v (x, x') J 1 ' (a/) , 



D+(S) 

dH r > 



TT c,I '{x')Gf If (x,x') - A P {xf){L v G H ) I '" (x,x') 



c'l' 



(2.12) 

(2.13) 
(2.14) 



We have changed the integration domain for A\ (x) from U = -D + (X) n D~(Y>') to D + (Y>) 
because Gf r (x, x') = unless x' G D~(Y!) by the assumption that x G U. The field Aj(x) 
given by Eq. ( 12.121) is the unique solution to the inhomogeneous equations ( 12.81) for the given 
initial data on E. 

In the next section we show how this formula can be used to reproduce the electric 
field in the future of the t = Cauchy surface in Minkowski spacetime with a static point 
charge. In Sees. HVj and [V] we show in examples that the retarded Green's functions for 
electromagnetism [26| and linearized gravity 17| in de Sitter spacetime generate the correct 
fields through this formula. 



III. ELECTROMAGNETIC FIELD IN MINKOWSKI SPACETIME 

In this section we demonstrate how the formula derived in the previous section works in 
the simple example of the static electric field compatible with a static electric charge in the 
future half of Minkowski spacetime. The Lagrangian density for the massless spin-1 vector 
field A a in the covariant gauge is given by 

C = -\F ab F ab - 1 {d a A a f , (3.1) 

where the electromagnetic field tensor is F ab = d a A b — d b A a . We adopt the Feynman gauge 
( = 1. As is well known, the retarded Green's function in this gauge is 

G a a ,(x, x>) = ^ ( Vf~'T X 'V - • (3-2) 

47T|X — x'l 

In Fig. [2] we show the support of the source field from a static point charge if the initial 
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Initial surface 




t = 



FIG. 2: A static charge q at the center producing the field in the shaded region 



surface is taken to be the t — hyper surf ace. Now, let a static point charge q be placed 
at the spatial origin. Then the corresponding 4-current is J a (x) = (g5 3 (x),0). With the 
retarded Green's function in Minkowski spacetime given by Eq. (13.21) we find the source field 
to be 



A a(S) { 



X 



4 f q s(t- a 



d x 



x'|) 



t'>0 

Q 

4/t|x 



47T 

■0(f-|x|)<5» 



x — X' 



6(t - t')5?5 3 (x' 



(3.3) 



where 9(t — t') is the Heaviside step function. We note that the Heaviside step function 
appears naturally due to the causal nature of the retarded Green's function. As illustrated 
in Fig. [31 the source field spreads as the time elapses. 

A* 



Time 




Fig. 3: The spreading of the source field as the time elapses. 
Next, we turn our attention to the initial field. We use the input field A a> = (A* , 0) with 



47TX' 



(3.4) 



where x' is the position vector on initial surface. This field is generated on the initial Cauchy 
surface by the charge before the time t' = 0. In the Feynman gauge ( = 1 we find that the 



Lagrangian density can be written as 

C= 1 -T abcd d a A b -d c A d , (3.5) 

where 

T abcd _ g bc g ad _ g ac g bd _ g ab g cd _ (3 

Hence the conjugate momentum current is 

vr ca = T cadb d d A b = -F ca , (3.7) 

where we have used the fact that the field (13.41) satisfies the Lorenz condition d a A a = 0. The 
only nonzero component of the past-directed unit normal to the hypersurface S is N t > = 1. 
We find from Eq. (12.141) that the initial field is given by 

A a(I) (x) = [ dX t , \&*'A* ■ G a a , - A a , (-d t 'G aa ' + d a 'G at ' - g t>a ' d y G ah '\] . (3.8) 
Jt'=o L V J\ 

Noting that A a > (—d t 'G aa ' + d a 'G at ') = and dropping the total spatial divergence 
El^diiA^G^), we find 

A^\x) = - [ MA'JLg** 

Jt'=0 en 

d*x> q W-^-^Sf 
47r|x'| 47r|x — x'| 

' nn *\-tWt- (3-9) 



47r|x 




Fig. 4: The behavior of the initial field as the time elapses. 

Fig. H] illustrates that the initial field eventually becomes zero for any given point as time 
elapses. From Eqs. (I3.3P and (I3.9P we find 

A a (x) = A a{s \x) + A a ^\x) 

' V, (3-10) 



47r|x 



recovering the correct field. The calculations in de Sitter spacetime are much more compli- 
cated, but we shall see that the initial field fills up the region that is not in the causal future 
of the source in a way similar to this example. 
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IV. ELECTROMAGNETIC FIELD IN DE SITTER SPACETIME 



The metric of de Sitter spacetime can be found by regarding this spacetime as the hy- 
persurface —(X ) 2 + ^ A i= i{X % ) 2 = H~ 2 in 5-dimensional Minkowski spacetime with the 
metric ds 2 M = —(dX ) 2 + ^ A i= i{dX 1 ) 2 . This hypersurface can be parametrized as X° = 
H~ l ta,nT, X 1 = H ~ 1 sec r cosx, X 2 = H~ x seer sin \ cos 9, X 3 = H^ 1 sec r sin \ sin 9 cos </>, 
X 4 = H' 1 sec t sin x sin 9 sin (J), with the range of the parameters given by r G (— tt/2, 7r/2), 
X, 9 G [0, tt] and <p G [0, 27r). The resulting metric, which is valid for the whole spacetime, is 



ds 4 



H 2 cos 2 r 



(-dr 2 + dx 2 + sin 2 xd& 2 ) 



(4.1) 



where 



dVl 2 = d9 2 + sin 2 9 d(f) 2 . (4.2) 

In this section we reproduce the electromagnetic field compatible with a charge q at x — 
(the "North Pole") and another charge — q at x = 71 (the "South Pole") using the retarded 
Green's function. [As is well known, the total charge on a compact space must vanish (see, 
e.g. Ref. [3l[). Hence it is not possible to have just one point charge in de Sitter spacetime.] 
The massless spin-1 Feynman propagator Q a a'{x,x') in the covariant Feynman gauge in 
de Sitter spacetime background, satisfying 

(-VfcV" + 3H 2 )Q aa ,(x,x') = -ig aa ,5\x,x') , (4.3) 

where H is the Hubble constant of de Sitter spacetime, has been calculated by Allen and 
Jacobson 26[. We need to give some definitions in Ref. 26] in order to present their result. 
Given a pair of spacelike separated points x and x' we define /i(x, x') to be the geodesic 
distance between them. The variable z is defined by 



= V a /i at x and n a 



^(1 + cos Hfi) . 



We also define the unit vectors n a 
the geodesic between x and x' . The vector n a is given by 

1 



(4.4) 

V a '/x at x' . These are tangent to 



n n 



2HJz(\-z\ 



d a cos H/j . 



(4.5) 



We also define the parallel propagator, g aa >, as follows: given a vector V a at x, the vector W a ' 
obtained by parallelly transporting V a to x' along the geodesic is given by W a = g a a V a . It 
can readily be seen that V a = g a a ,W a . The parallel propagator is given explicitly as 



9a 



77j ( d a d a r cos H[i d a cos Hfi ■ d a > cos H[i 



H 2 



2z 



Then, the Feynman propagator of Allen and Jacobson is 

Qaa'(x,x') = a(z)g aa >(x,x') +/3(z)n a (x) n a >(x') 

where 



a(z) 
(3{z) 



H 2 
48^ 

H 2 
24T 2 



3 1 
+ - + 

1 — z z 



2 1 
- + -s 



log(l - z) 



1- - + 

z 



1 1 . 

- - - log(l - z) 



(4.6) 
(4.7) 

(4.8) 
(4.9) 
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The Feynman propagator Qaa>(x,x') for points that are not spacelike separated is defined 
by analytic continuation. It is analytically continued around the singularity at /i 2 = 0, i.e. 
z = 1, by replacing [i 2 by fx 2 + ie, where e is an "infinitesimal" positive number. As a result 
1 — z is replaced by 1 — z + ie. 

The retarded Green's function is proportional to the difference of the values of the two- 
point function across the branch cut [1, oo) on the real axis. (The relation between various 
Green's functions can be found, e.g. in Ref. [32| for scalar fields. The Green's functions for 
bosonic higher-spin fields are related in the same way.) We note that the singular terms at 
z = in (3{z) cancel out. The retarded Green's function G aa >(x, x') is given by 

G aa > (x, x') = -i9(t - t') [Q aa > (z + ie) - Q aa ' (z - ie)} , (4. 10) 

where r and r' are the time coordinates of points x and x', respectively. The explicit form 
of G aa '(x,x') can be found by using the Sokhotski formula, 

' ~ 1 =p 17T<J(1 - z) , (4.11) 



1 — z ±ie 1 — z 
where P denotes the principal part. We also use 

log(l - z ± ie) = log |1 - z\ ± m9{z - 1) . (4.12) 

The resulting retarded Green's function is 

G aa ,(x,x') = G^,(x,x')5(l -z) + G^(x,x')9(z - 1) , (4.13) 



where 



= ^9aa>(x,x'), (4.14) 



i x i x ) — 



12tt 



z + 2^2 ) gaa '( x,x '} + ~2~ n a(x)n a >(x') 



(4.15) 



Here, the factor 9(r — r') is understood though not explicitly written since it will not affect 
our calculations. Note that the singularity of n a n a i as z — > 1 [see Eq. (14.51) ] is canceled by 
the factor z — 1. 

A solution to Maxwell's equations with the two charges at the Poles is found most readily 
in the static coordinate system with the metric 

ds 2 = -(1 - H 2 R 2 )dT 2 + (1 - H 2 R 2 )' l dR 2 + R 2 dtt 2 . (4.16) 

This coordinate system is related to the global coordinate system (14. ip as 

1 , cosx + sinr 
in cos x — sin r 

R = -P*-. (4.18) 

H COST 



The static electric field with charge q at the origin R = is given by 

rpTR = g 

4ttR 2 ' 



(4.19) 
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This equation is solved by 

with all other components vanishing. We have added a constant qH/4n to a more natural 
solution At = —q/AirR to make the scalar quantity A a A a non-singular at the horizon R = 
1/H. This field satisfies the field equation (□ — 3H 2 )A a = at R 7^ and the Lorenz gauge 
condition V a A a = 0. Hence it should be reproduced by the retarded Green's function (14.131) . 
In the global coordinate system this field is 

q cos r cot v 

A T = -f 4.21 

47T cos t + sin x 

A x = -f SinT . (4.22) 
47T cos r + sin x 

Under the transformation % 1— > 7r — which interchanges the North and South Poles, the 
field A a simply changes its sign. Thus, it is clear that there are two charges, q at the North 
Pole and — q at the South Pole. 

Our calculations to reproduce the field (14.201) by the retarded Green's function can be 
performed most simply in the conformally-flat coordinate system with the metric 

ds 2 = -^(-dX 2 + dr 2 + r 2 n 2 ) . (4.23) 

This coordinate system covers only the upper half of de Sitter spacetime, which is the causal 
future of the charge at the North Pole, represented by the shaded region of the Carter- 
Penrose diagram in Fig. [TJ This coordinate system is related to the static coordinate system 
by 

T = -^\og[H 2 (X 2 -r 2 )} , (4.24) 

Note that the time variable A increases toward the past. The nonzero components of the 
electromagnetic potential in the new coordinate system are 

* = sf(?-tta)' (426) 

^ = "S7Ta- (427) 

We need to express cosH[i(x,x') in terms of the coordinates of the two points x and x' in 
order to find the explicit expressions of 2, n a and g aa > through Eqs. (14.41) . (14.51) and (14.61) . 
By writing down H~ 2 cos Hfx = — X°X ' + Ylt=i X % X 1 ' in conformally-flat coordinates we 
find 

A 2 + A' 2 - r 2 - r' 2 + 2rr'[cos 9 cos 9' + sin 9 sin 9' cos(0 - d>')} 
cos Hfi = L__ K^ni . (4.28) 

The function cos Hfi is extended naturally to the case where the two points are not spacelike 
separated by adopting this formula as it is. 
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Now, we calculate the source field from the charge q at the North Pole. We adopt the 
convention that the unprimed coordinates are those of the point where the field is calculated 
whereas the primed coordinates are those of the charge. The 4-current corresponding to the 
charge at the North Pole can readily be determined by noting that it must be a conserved 
current. It is given by J a = ( J x , 0), where 

J A ' = -g# 4 A' 4 5 3 (x'). (4.29) 

Here the delta function is defined to satisfy 

dvyy^x') = 1, (4.30) 



where dVVT 7 = dr'dO'dft r' 2 sin 6'. We have from Eqs. (14.281) and f)4.29p 

^g'5{\ - z)J x ' = -?2^-5(\' - A - r) v9<*V) > ( 4 -31) 



y^e(l-z)J x ' = -q9(X'- A-rjVy^x'). (4.32) 

We first cut off the current at A' = Ao by using the 4-current J a ' 9{Xq — A') and then take 
the limit Ao — > oo in the end. (Recall that A' increases toward the past.) Substituting the 
retarded Green's function (I4.13P in Eq. (12.131) and using Eqs. (I4.3ip and (14.321) . we find 



OO 

dX G aX , 

A'=A+r -'A+r 



9(X - A - r) , (4.33) 

=o 



4 S, = ^fz-r 1 T + ^)^.-A-r), 



and hence 



An \r r + A 3A 

with Ag S ^ = = 0. The source field vanishes for r > A — A, i.e. outside the causal 
future of the world line of the truncated charge J a '9(X — A'), but in the limit A — > oo the 
field becomes nonzero everywhere in the causal future of the charge at the North Pole, which 
coincides with the part of de Sitter spacetime covered by this conformally-flat coordinate 
system. Thus, the Heaviside function can be dropped in this limit. We note that for 
r < Ao — A the field A a does not agree with that given by Eqs. (I4.26P and (I4.27P nor does 
it satisfy the Lorenz gauge condition because of the term 1/3A in Eq. (14.341) . Nevertheless, 
the source field gives the correct field strength F a b because the extra term is of the form 
V a A with A oc log A. We will show that this extra term will be canceled by the initial field. 

Next we turn our attention to the initial field. The coordinate system we have used, 
which covers the causal future of the charge q at the North Pole, is not appropriate because 
past infinity, where we need to specify the initial data, is not covered at all. For this reason, 
we will work in the coordinate system which covers the causal past of the charge. Thus, we 
consider the coordinate system covering the region illustrated in Fig. [5j for which the metric 
takes the same form as before: 

ds 2 = —L- i-dX 2 + df 2 + f 2 dQ 2 ) . (4.36) 
H 2 X 2 V / 
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Boundaries of coordinate patch (A = oo 




increases 



Surfaces of constant time A 



r increases 



* — the world line of the North Pole 



FIG. 5: Coordinates (A, f, 6, 4>) covering the lower half of the manifold. 



The coordinates A and f are related to A and r for the upper-half spatially-flat coordinate 
system by 



A 



A 



H 2 (\ 2 -r 2 ) 

f 

H 2 (\ 2 -r 2 ) 



(4.37) 
(4.38) 



in the overlapping region. The angular variables are related trivially as 9 = 9 and <fi = <f). 
The time coordinate A 6 (0, oo) runs forward from the infinite past (A = 0) to the infinite 
future (A = oo). The initial surface will be a hypersurface of constant A , which will be left 
arbitrary though we are mainly interested in the limit A' — > 0. 

We let the unprimed coordinates correspond to the point where the field is calculated and 
the primed coordinates correspond to the point on the initial surface S. In this coordinate 
system we have A > A'. We find from Eq. (14.201) that the input field on the initial surface is 
given by 



4tt \f' f' + A' 



(I 



4vr f + A' 

with Ax, = At, = 0. The surface element is 



(4.39) 
(4.40) 



dY,N c , 



f' 2 sin 9' 



H 4 \ 



dr'd9'dd)' 8H . 



(4.41) 
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We only need the A'-component of the conjugate momentum current, ir x ' a ' . In the Feynman 
gauge ( = 1 the only nonzero component of the conjugate momentum current is 



q H A \' 4 
47r r' 2 



n X _ g\ A g r f {d , x/Afi _ dfjAxi) = -j-^--. ( |. |2 ) 

The differential operator acts on the Green's function as 

{L n G) a c ' a ' = -±-{H 2 [5'(z-l) + 25(z-l)}g c ' a, d a cosHfi 
lb7r i. 

+25\z - l)d a d [c ' cosHfi -d a ' ] cos . (4.43) 

The function cos Hfi is given by Eq. (I4.28P with all variables replaced by the hatted equiva- 
lent. (For example, A is replaced by A.) The following formulas from Ref. [26[ are useful in 
deriving Eq. (14. 43 ft : 

V B 'Z = -Hy/z(l-z) n a ,, (4.44) 
V a >n b > = H< f Z l \ {g a >b< - n a >n v ) , (4.45) 

2y/z(l ~ Z) 

V a /n a = (g aa , +n a n a >) , (4.46) 

2s/z{l - z) 

^a'9ab' = H^^—^(g a , b ,n a + g a > a n b ,) . (4.47) 

We specialize to the point (A, f, 8 — 0) — <j> is undetermined — for calculating the field 
without loss of generality because of the spherical symmetry. The integral on the initial 
hypersurface with the time parameter A' in Eq. (12.141) . which gives the initial field, is over 
f and 8'. (The integration over 0' gives a factor of 2tt.) Due to the causal nature of the 
retarded Green's function the nonzero contribution comes only from the intersection of £ 
with the causal past of the point (A, r, 8 = 0), which is a 3-dimensional ball, B. There 
are two cases to be considered separately: i) f > A — A' [(A', 0) ^ B] and ii) f < A — A' 
[(A , 0) GB]. For the case f > A — A' (see Fig. [6]), for given 8' the range of integration for 



the radius f' is from r /_ to f' + , where f /=t = r cosd' ± y (A — A') 2 — f 2 sin 2 8'. The variable 8' 

is integrated from to = arcsin ■ For the case f < A — A' (see Fig. [7]), the integral 

over f' is from to r' + for fixed 8' . Then the angle 8' is integrated from to n. Using 
Eqs. ( 14TT4D . (JHSD, (14^391) . pg| . (Q2|) and (031 in Eq. (12H|) . we obtain the initial field 
in the part of the spacetime covered by the lower-half spatially-flat coordinate system as 



A f = r 1 i x.. vi 

4tt r + A, 



5 



VA A + f A-f/ 



12tt VA 
Q 

An(f + X)' y ' ' 12tt VA + f A 



0(A — A' — f ) , (4.48) 



4° = ^— -^(^ - A + A') + -f - f — !— - — ) 8{X - A' - f) , (4.49) 
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Field point 
(A,M = 0) 



spatial origin 
x' = 



world line of a 
' static charge 




(A - A') 2 - r 2 sin 2 §' 



FIG. 6: The causal past of the point (A, f, 6 = 0) in the case f > A — A . 



Field point 
(A,M = 0) 



A -A' 




spatial origin 
x' = 



world line of a 
static charge 



(A - A') 2 - r 2 sin 2 B> 
FIG. 7: The causal past of the point (A, f, = 0) in the case f < A — A'. 

with = = 0. The initial field in the reerion with f > A — A' is equal to the field 
given by Eqs. (14.391) and (14.401) as expected because this region is not influenced by either 
charge. In the limit A' — > the terms proportional to 9{\ — X' — r) gives the field A ( P m 
the region covered by both (A, f) and (A, r) [see Eqs. (I4.37P and (14.38P ]. This field can be 



15 



expressed in the upper-half coordinates (A, r) as 



overlap 12ttA v ' 

with all other components vanishing. This term cancels the extra term in Eq. (I4.34p . making 
the sum agree with the expected field in this region as well. 

The calculation in this section shows that the retarded Green's function for the massless 
vector field obtained from Ref. 26[ correctly reproduces the field compatible with a charge 



q at the North Pole and a charge — q at the South Pole in the causal past of the charge q 
at the North Pole. It is clear that this is also the case in the causal past of the charge — q 
at the South Pole. Then, the uniqueness of the solution to the gauge-fixed field equations 
with given initial data on a r-constant surface with r < in the global coordinate system, 
which is in the causal past of the two charges, implies that the retarded Green's function 
reproduces the field correctly over the whole spacetime through Eq. f)2.12p . 

V. GRAVITATIONAL FIELD IN DE SITTER SPACETIME 

In this section we reproduce the linearized gravitational field compatible with two point 
masses M, one at the North Pole and the other at the South Pole of de Sitter spacetime, 



using the retarded Green's function obtained from the two-point function found in Ref. [17 
(One cannot have a single mass point in de Sitter spacetime because the total conserved 
charge corresponding to a de Sitter boost symmetry must vanish. This fact is related to the 
linearization instabilities of the Einstein equations about any spatially compact spacetime 



with Killing vectors such as de Sitter spacetime [3ll. l33l|.) 



A. Linearized gravity for a point mass in a covariant gauge 

We first construct the field to be reproduced starting from the Schwarzschild-de Sitter 
solution in the static coordinate system: 

ds 2 = - (l - H 2 R 2 - ^ dT 2 +(l- H 2 R 2 - ^ 1 dR 2 + R 2 dn 2 . (5.1) 

As is well known, this solution represents two black holes of equal masses M in de Sitter 
spacetime. If we write this metric as g^ u = g a b + h a b, where g a b is the background de Sitter 
metric (I4.16p . then the linear term h a b is singular at the horizon R = H^ 1 . However, after 
the coordinate transformation given by 

R - U-ImH^R-^, (5.2) 
T -> (1 + 2MH)T, (5.3) 



16 



we find 



hij 



2M 



(1 - HR) (l-HR- ^ h2 R 2 ^J 



2M 



3R(1 - H 2 R?) 



1 + HR 



-MR{1 + 2HR) Vij , 



(5.4) 
(5.5) 
(5.6) 



where i, j = 6, <p and where rjij is the standard metric on the unit 2-sphere. The singularity 
of the component at the horizon R = H^ 1 is a coordinate singularity, and this field 
is non-singular at the horizon as one can see by expressing in the global coordinate 
system. The expression of h a b in global coordinates also makes it clear that there are two 
mass points, one at the North Pole and the other at the South Pole. 
The trace and divergence are 



h a 



Hence h ab satisfies the gauge condition 



AM 

~~3R J 
10M 

3R 2 a 



(5.7) 
(5.8) 

(5.9) 



with (3 — |. This gauge condition is exactly the one for which an explicit form of the 
Feynman propagator is given in Ref. (l7| . 

It is possible to make the field traceless and divergence-free at R ^ by a further 
gauge transformation, 

R ^ R -^r^~ H2r2 )- ( 5 - 10 ) 

However, the resulting i?i?-component will behave like R~ 3 at R = 0. We do not adopt this 
solution for this reason and will work with that found above. 

In the conformally-flat coordinates covering the upper half of the spacetime we find 



fl/i A 



AM 
3H 
AM 
3H 
AM 
3H 
AM 



1 



1 



2rA 



r(r + A) (r + A) 2 

2 1 2 

r(r + A) (r + A) 2 A 2 

2 1 2 



2rA 



r(r + A) (r + A) 2 

^,2 ■ 



rA 



3H V 2A 



+ p J Vij, hJ = o< 



(5.11) 
(5.12) 
(5.13) 
(5.14) 



with all other components vanishing. We reproduce this field using the retarded Green's 
function obtained from the result of Ref. [l7| with /3 — |. 
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B. The retarded Green's function for linearized gravity 

The Lagrangian density of pure gravity with positive cosmological constant 3H 2 is given 

by 

£ fuI1 = (167T)- 1 y/=§(R - QH 2 ) , (5.15) 

where g is the full metric and R is the Ricci scalar. We write the metric as g ab = g ab + h ab , 
where g a b is the background de Sitter metric, expand the Lagrangian (15.151) to second order 
in h ab , and then include the gauge fixing term, 



C 



gf 



-g 



32ira 



ab l+ P V * h 



V a h a0 - 







(5.16) 



where a and (3 are the gauge parameters. Specializing to the case with (3 = |, we find the 
Euler-Lagrange field equations to be 



— L a b cd hcd 



ph ab - 

25 
4a 



1 

2a 



+ 
0. 



{V a V c h c b + V b V c h c a ) + 

5' 



5 

2a 



V a V h h 



a 



g ab Dh + \g ab 1 - - V c V d h cd - H 2 {h ab + \g ab h) 



(5.17) 

(The definition of L ab cd differs from that in Ref. 17] by a minus sign.) The Green's function of 
the operator L ab cd given by Eq. (15.171) was calculated in Ref. [l7J by analytically continuing 
the corresponding Green's function on S 4 of radius H~ l . The corresponding Feynman 
propagator Q aba ' b '(x, x') satisfies 



Lt cd Q cda < b ,{x,x' 



-i5 



ab 



a'b' 



(x, x) 



(5.18) 



where the subscript x in indicates that the differential operator L . acts at x rather 
than x'. Here the 5-f unction is defined by 

^g~W) d 4 x' 5 aba , b ,(x,x')f a ' b ' \x') = f ab (x) (5.19) 

for any smooth compactly-supported symmetric tensor f ab . As in the original graviton two- 
point function in the gauge a = 1 and f3 = —6 16j, the covariant graviton propagator in the 
gauge (3 = | can be expressed in terms of bi-tensors formed using the five basic (bi-)tensors 
gab,9a'b',g a a',n a ,n a >. They are 



^aba'b' 

0(2) 
w aba'b' 

0(3) 

o 
o 



(5.20) 
(5.21) 
(5.22) 
(5.23) 
(5.24) 

where g ab > and n a are defined by Eqs. (14.51) and (I4.6p . and where g ab and g a > b / are the metric 
tensors at x and x', respectively. The covariant graviton propagator is 



(4) 

aba'b' 
(5) 

aba'b' 



gabga'b' , 

9aa'9bb' + 9ab'9ba' , 

gabn a >n b > + g a ' b 'n a n b , 

9aa>n b n b > + g ab 'n b n a , + g bo! n a n v + g bb >n a n a , 
n a n b n a 'n bl , 



Qaba'b'(x,x') 



H 2 

167T 



(5.25) 



k=i 
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where f(k)(z) are scalar functions of z(x,x') given by 



/(!)< 



22a 



15 



4 
9 



17a \ 1 



45 



a 



( 45 9 z 2 



f(2)(z) 



+ 

2 
3 

+ 

4 
9 



4a 4 



a 



14 



a 



4 



13a 



2a 



1\ 1 



45 



a 



45 9/z 



1 

3 
a 
45 
5a 
9" 



4a \ 



4 



log(l - z) , (5.26) 



3/1-2 



6a 



a 



2a 

T 



/(4)(*) 



/(5)0 



+ (!-«) 

2 16a 

3 ~ ~1T~ 

+ (!-«) 
8 152a 



9 
31a 
"45~ ~ 
a 1 
5 ~ 3 
6a 
"5 
3a 

T 
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2- — + --- -T+ tt-t: 



a 



4 



+ 



4 



+ 



-2 • 



2 

7a 
45 
2 

2 



/7a 
V 45 



45 9 
4 a\ 4 
3 ~ 15/ I2 
2 



2 J 
+ 



log(l - z) 
la 1 



(5.27) 



9 1-2 



+ 



a. 



4 



9 

+a 



15 



24a 

+ — z + 

1 a 
3 ~ 5 



4 /4 

a 1\ 8 
15 ~ 9J 1? 
4a\ 8 



a 

45 ) ~z 2 



15 y 2 3 

a — 33 



log(l - z) 
1 



(5.28) 



9 



1 



+ 



- - + 



+ 



a 

15 / 2 ' V45 
a _ 4\ 16 
45 ~ 9/ z 1 



4 a 
9 ~ 45 J : :! 
4\ 16 



loefl 



1 a 
0/ : 2 V3 + 9 



- z 
32 



(5.29) 



log(l 



(5.30) 



In finding the retarded Green's function it is important to note that the products n a rifc and 
n a n a i contribute a factor (1 — z)~ l [see Eq. ( 14.51) ]. which is singular as 2 — > 1. (We did 
not face this issue in the electromagnetic case because the bi-vector n a n a i was multiplied by 
1 — z in the two-point function.) We note that there is no term of the singularity structure 
of the form (1 — z)~ n \og(l — 2), n > 1, in Q a ba'b'(x,x'). Such a term would have caused 
difficulties in finding the retarded Green's function. 
It is convenient to define 



C (i) 



abc'd' 



9ab 9c' d! 



(5.31) 



C {2) abc'd' = jp {9ab d c > cos H fi ■ d d > cos H fi + g dd , d a cos Hfi ■ d b cos Hfi) , (5.32) 



C (3) 



abc'd' 



H 4 



(d a cos Hfi ■ db cos Hfi ■ d c i cos Hfi ■ 841 cos Hfi) 



(5.33) 



The bi-tensors C^ l J bd , 1 n = 1,2,3, are defined by swapping b and d in these formulas. 
The retarded Green's function G a b a 'b'{x,x') can be obtained in the same manner as in the 
electromagnetic case by noting 

G aba , b ,(x,x') = -%9{t - t') [Qaba'b'(z + i"e) - Qaba'b'{z - ie)] , (5.34) 

and using the Sokhotski formula (14. lip and its derivatives, 



1 — z ± ie 



n+l 



1 

n! 



d n 
dz n 



±(-\) n+1 m5 {n \l- z) 



(5.35) 
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and Eq. (I4.12p . We obtain the retarded Green's function in this manner as 

H 2 



where 



G 



(o) 

aba 1 1 



^ aba'b' 



G 
G 



(2) 

aba'b' 
(0) 

aba'b' 



G 



aba'b' 



in 



G { Z V 6(1 -z) + G& v 6'(l -z) + Gfl, b ,5"(l - z) 

+G { 2 a , b ,e{z - 1)" , 



(5.36) 



— ^ C^aba'b' + (— T^j (C^ aa 'bb' + C^aWa) + g 1^ 'aba'b' 

-1(19 + 13a) {CW ab , ba , + C^ aalw ) - ~ (l + ^ C^ aha , b , , (5.37) 



- + — ) C^ aba ' b r 



(2z-l) 
3a 



3a \ 2 /4 14a \ 1 /2a 8\ 1 
~T/3z + V3~^5/^ + V45 ~ 9 ) i 3 



aba'b' 



+ 
+ 
+ 
+ 



5 X J l 2+ 9z 2 ) + V9 45) z 3 



3a 

T 

3a 



' ~ ~5~ J 3^2 + I 3 ~ ~45~ / z 3 + V 15 3 / z 4 



aa'bb' 
4\ I 



1 2 \ /2a 8\ 1 
2z 2 + 9z^ + V 45 ~9 i* 



r<( 2 ) 

ab'a'b + C {2) a a'bb') 



3a \ 1 / 4 _ a 
T ~ /6^ + V9~45/;z 5 



(5.38) 
(5.39) 



(5.40) 



C. The source field 



The formula (12.121) for linearized gravity is 

h ab (x) = h a S b \x)+h < f b \x) 



where 



,T / d 4 x'^Gaba'b'(x,x')T a ' b '(x') 
'D+(S) 



n c ' a ' b ' (x')Gaba'b' (x, x') - h a ' V {x') {L^G) ab c a h (x, x') 



(5.41) 

(5.42) 
(5.43) 



We explicitly show that this formula reproduces the field h ab in Sec. IV Al compatible with 
two mass points M at the North and South Poles. In the conformally-flat coordinate system 
covering the upper half of the spacetime, the corresponding stress-energy tensor is propor- 
tional to 5 3 (x), representing the mass point at the North Pole. It can readily be found by 
using the conservation law V a T ab = and dimensional analysis as 



T AA = MH b X b 5 :i (x.) 



(5.44) 
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with all other components vanishing. 

The source field (15.421) is found using Eqs. (I5.36P and (I5.44p as 



MH 



dX' X' 



G% f>f 6(l-z)+G&,x#{l-z) 



(6) 



(5.45) 



Again we cut off the stress-energy tensor at A' = Ao by multiplying it by 6*(Ao — A') and let 
A —>■ oo in the end. By integrating by parts to reduce the order of derivative of the delta 
functions and letting J Q A ° dX'6(z — 1) = J x °dX', we obtain the AA-component as 



h 



(S) 
AA 



MH J 



dX' 



A AAA' A' + 



d / d 



TTT7 I A'G?> I ^ ] I [ ^ 



dX' \ dX' 

Ao 



'AAA' A 



d_ 

dX> 



dz 



A AAA' A' 1 Q^, J 



dz\ 1 



-r 



f dz 



+ 



dX' X'G 



AAA' A' 



AM 
3H 



X+r 

2 



dX' 



1 9a 
+ 



\dX f 



-r 



6(1 -z) 



_r(r + A) (r + A) 2 2rA 5A 2 



9(X -X-r). 



(5.46) 



The other components can be found in a similar manner as 



i {S) 

VA 



4M 
3H 
AM 



AM 
3H 



r 

2A 



1 



1 


1 


(r + A) 2 


2rA 


1 


2 


(r + A) 2 


rX 


3a\ r 2 ' 





i 3a A 1 

|2 + T>- 



9ar 



A 2 

e(X -X-r) 



A 2 



^^(Ao - A - r) 



fl(Ao-A-r), (5.47) 
(5.48) 
(5.49) 



( s) 

with the remaining components vanishing. The source field h ab agrees with the field we 
found earlier given by Eqs. f l5.11l) -( l5"TT4l) if we let a = 0. We also find that the a-dependent 
contribution can be expressed as V a A 6 + V&A a where 



QaMHX 



5 

2aMHr 



hgHX 



l-3\ogHX) 



(5.50) 
(5.51) 



with A e = A* = 0. Thus, the retarded Green's function reproduces the field correctly in the 
causal future of the point mass up to a linear gauge transformation. The source field with 
q^O does not satisfy the gauge condition: 



F a = V b h ba - -V a h 



-6aMH 3 X 
-QaMH 3 r 




a = X . 
a = r , 
,a = 9, 



(5.52) 
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The vector F a is a Killing vector, and as a result the contribution of the gauge-fixing term to 
the field equations, a^ 1 (V a Fb + V bF a — bgabV C F C ) , vanishes. Since this could not be the case 
if the spacetime had no Killing vector, we expect that the source field by itself would not 
satisfy the linearized Einstein equations for a ^ in a spacetime without a Killing vector. 

In the next section we will show, with a = |, that the initial field given by Eq. (15.431) 
reproduces the field in Eqs. (I5.11l) -( l5.14l) correctly at points which are not in the causal 
future of the mass points and exactly cancels out the extra pure-gauge term from the source 
field in their causal future. 



D. The initial field 

The calculation of the initial field is rather complicated, and we specialize to the gauge a = 
|, which somewhat simplifies the retarded Green's function. Again we use the conformally- 

flat coordinate system (A, f, 9, 0) covering the causal past of the point mass at the North Pole. 
We let the input field h a ty on the initial hypersurface A = A be equal to Eq. (I5.11l) - (l5.14p 
with each variable changed to its hatted equivalent. 

Let us first find the conjugate momentum current ix cah . As in the electromagnetic case 
we denote the coordinates at the initial hypersurface by primed variables, (A', r', 6', 0'), and 
those at the point where the field is calculated by unprimed variables, (A,f,#,0). The 
Lagrangian density for linearized gravity with gauge-fixing term with a = | and (3 = | can 
be written as 



C = [T abcd£ fV a h bc V d h ef + S abcd H 2 h ab h cd ] . (5.53) 



where 



T abcdef = ^ g ab g cf g de _ }_ g ad g be g cf _ ]^ g ad g bc g ef + g *c g bd g ef + g df g ae g bc ? ^ ^ 
gabcd = _ g ac g bd + igOb g cd _ (5 55) 

The conjugate momentum current is given by 

1 dC 



„cab 
71 



V~9 diy c h ah ) 

= 1 ^cabdef + T cbadef^ 

= -\V c h ab -\g ab V c h + 2g c{a V b) h, (5.56) 
where we have used the covariant gauge condition to replace Vbh b a by \V a h. As for the 
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conjugate momentum current we only need the components Tc x ' a ' b ' . They are 

AMH 5 



71 



A'A'A' 



7T 



X'f' 



X'X'f' 



3 

AMH 5 
3 

AMH 5 



-A 



-A 



r-i 



A 



A' 2 




1 


1 


(A' + f') 3 


A' 


+ r' 




A' 2 




1 




(A' + r') 3 


A' 


+ r' 


• 


A' 2 


- + - 


1 


A' 



(A' + f') 3 A' + f 



7T 



3f' 3 



Next, we evaluate 



c'a'V 



ab 



1 I ' rpc'a'b'd'e' f , r r c'Va'd'e'f'\ V7 

5# Vd'Crafc — ^ V Cr a ft + <? 
+g Vd'^ab — —g V (j a b , 



ab 



(5.57) 
(5.58) 
(5.59) 
(5.60) 



(5.61) 



where G a b = G a b a Recall we have expanded G a bc'd' in terms of the tensors C^ c , d , and 
Cac'bd'i i = 1)2)3, defined by Eqs. (I5.3ip - fl5.33p . We need the covariant derivatives of these 



tensors in order to find (L 7r G) ah cab . These are evaluated in Appendix IB1 

As in the electromagnetic case we group together the terms proportional to 5(1 — z) 
5'(1 - z), 5"(1 - z), 5"'{l - z) and 6(z - 1) in (L W G) 



c'a'V 



{LttG) 



c'a'V 



ab 



IP 



s. 



(or 



5(1 -z) + " ~ S'(l -z) + S%T - " 5"(1 - z) 



ab 
'a'V 



as 



ab 



(2) c 



+S^-5">(l-z) + S^ a r D 6(z-l) 



c'a'b 



c'a'V 



(5.62) 
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where 



c'a'b' 



ab 



c'a'b' 



ab 



and 



c'a'b' 



S { 9 ) 



4 /~~Ti \ I 1 i c'a' b' 25 „/ a! V H a'V d da 1 V 

- Vz(l - z) I llg c g [a b n b) - —g [a c g b) a n b - —g a g (a c n b) - g ab g c n b 



11 _/ L/ J 41 _/ , . 

+ Y9(a 9b) n c + — g a g ab n c -z{l-z) 



165 „,„ 



112 , , 



g a n a n b n c + — # c a n a n b ri 



b' 



130 / / u 152 / i/ / 11 / 
+— 9(a n b) n a n° + —g {a n b) n° n c + —g a bn a n 



a'^b' n c' 



(5.63) 



-274z 2 (l- z) 2 n a n b n a 'n b 'n c '\ , 



4 
9 

11 n' V r-' 35 „/!/ „/ , . 



247 -i 41 /_/ 1 

— g a n a n b n c + — c/ c n a n b n b 
6 3 



+27g {a c 'n b) n a 'n b ' + 3Qg( a a ' n b) n b ' n c ' + 202^(1 - z)n a n b n 



a 'n b 'n c ' 



(5.64) 



±y/z(l-z)z(l - z) I -^/<>'n n 6 n c ' - 28 (/"'n^n 6 ' + ^ (a c 'n fe) n a 



' 6' 



-A7g {a a 'n b) n b 'n c ' + hg ab n a ' n b V - 184z(l - z)n a n b n a V V } , (5.65) 

(5.66) 
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z 2 (l - z) 2 y/z(l — z)n a n b n a ' n y n°' 



c'a'b' 22y/z(l - 

ab ^ ~ A 



9z 4 



-2 I 1 + y ) /°'<? ( /n 6) - 2< 7(a c V'n 6 ' - g { /g b) b 'n c ' 



+ ('f - - l) ^/<V' + (f - l) / - (x - 2 " + [) ^ W 



+2(1-2) 



(\ - # a ' b 'n a n fe n c ' + (2 - z)g ab n a ' n b ' rf' + 2g {a c 'n b) 

2„ „ „a' 6' 



+ 4# (a a 'n b) ri 6 V + 2 ( 1 + y ) g c ' a ' n a n b n b ' 



2z 



- 4(1 - zyn a n b n a n b n c } .(5.67) 



These formulas have been simplified by using the identities 

'(-l) fe+1 n! 



[/(*)- /(I)] *<»>(!-*) = £ 

fc=i 

n r 



fc!(n-Jfe)! 



fc=0 



ri! 



fc!(n - k)\ 



f k (z)5 {n - k \l- z) 



f {k \l)S (n ' k \l- z) 



(5.68) 
(5.69) 



To evaluate the initial field we substitute in Eq. (I5.43P the initial data h ab and -n x ' ab 
given by Eqs. (15. lin - (15 . 14p (with hatted variables) and Eqs. fl5.5Tn - f15.60n . respectively, and 
the Green's function G ab a' b ' given by Eqs. (I5.36f) and the function (L n G) ab a given by 
Eq. (15.621) . In performing the integral over the initial hypersurface it is easier to integrate 



24 



over the angle 6' first and then over p. For this purpose we re-express the delta function as 
f(1 -, ) = ^,f aM »-''+^-*) , V (5.70) 



2f f ' 



This can be used to perform the ^'-integral of the form 

dp f d cos 9' f(P, cos §')6 {n) (l - z) 
o J-i 

for n — 0, 1, 2, 3, where /(f cos #') is any function of P and cos 9'. The outline of the lengthy 
calculation together with some intermediate formulas used are given in Appendix [Bl The 
resulting initial field is given by 



h 



[i) 
"XX 



AM 
3H 

AM 



f(f + A) (r + A) 2 2fA 
2 1 



h 



(I) 



3H 

AM 
3H 

AM 



+ 



9{f - A + A') 
2 1 



3 

+ — 



r(f + A) (r + A) 2 f(f - A) (f - A) 2 A 2 
2 12 1 



0(A- A' -r), (5.71) 



2 

+ — • 

f(f + A) (f + A) 2 A 2 2fA 



1 



6{r - A + A') 



Af 



3i7 

4M 
3¥ 

AM 



1 



f(f + A) (f + A) 2 r(f — A) (f — A) 2 
2 1 2 



5 

I 2 



0(A - A' - r) , (5.72) 



f(f + A) (f + A) 2 fA 
2 1 



(f - A + A') 



3# 



+ 



1 



3f 



r(f + A) (f + A) 2 f(f-A) (f-A) 2 fA 2A 2 



h 



x0(A-A'-f), 
4M 



2 1 



r r 
2A I 2 



0(f - A + A') - 



4M f 2 



0(A — A' — f ) Vi 



(5.73) 
(5.74) 



We see that the input field on the initial hypersurface with the time variable A' is reproduced 
for f > A — A', i.e. in the region which is not in the causal future of either of the two mass 
points. 

As we observed before, the initial Cauchy surface coincides with past infinity in the limit 
A' — > 0. In this limit the terms proportional to 6(X — A' — f ) = 6(X — f ) give the initial field 
in the overlap region covered by both the upper- and lower-half conformally-flat coordinate 
systems. By transforming the initial field in this region from the lower-half conformally-flat 
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coordinates system to the upper-half one by Eqs. (14.371) and (I4.38p . we find 



h (I) 
xx 

him 



a Xr 



overlap 



overlap 



overlap 



overlap 



4M 3 


(5.75) 


6H A z 


AM 1 


(5.76) 


3H A 2 ' 


AM 3r 


(5.77) 


3H 2A 3 ' 


AM r 2 


(5.78) 


3H X^ J ' 



where i,j = 8,<f). By comparing these formulas with the source field given by Eq. (I5.1ip - 
(I5.14p and recalling that the initial field is computed with a = I, we find that the extra 



contribution in the source field h ab is canceled exactly by the initial field h^ b ' . Thus, the 



3 ' 



(J) 



field h r , 



h[ S J + h„i agrees with the input field for a = | in the causal past of the North 



Pole. The calculation of the field in the causal past of the South Pole is exactly the same. 
Hence, by the same argument as in the electromagnetic case, we see that the field h a b is 
correctly reproduced by the retarded Green's function through the formula (15.41 j) . 



VI. SUMMARY 

In this paper we showed in examples that the retarded Green's function does repro- 
duce the electromagnetic and gravitational fields in de Sitter spacetime through the formula 
(I2.12p . This verification is significant because it illustrates the fact that the de Sitter in- 
variant construction of the retarded Green's function is applicable to generate the fields 
which satisfy the equations of motion. Thus, the retarded Green's function should not be 
abandoned contrary to recent claims. Our calculations also serve as a consistency check for 
the propagators for electromagnetism [26] and linearized gravity [TtJ in de Sitter spacetime. 

It is interesting that the time component of the electromagnetic field coming from the 
charge at the North Pole in de Sitter spacetime possesses an extra term 1/3A which does not 
satisfy the Lorenz gauge condition. This term is canceled out by the initial-data contribution 
in the overlapping region of the spacetime (and thus the Lorenz gauge condition is satisfied). 
This result illustrates the need for including the initial data on past infinity satisfying the 
Gauss constraint equation in a spacetime with spacelike past infinity. (The field at a point 
which is not in the causal future of either charge is reproduced exactly by the field from the 
initial data alone.) 

The way the linearized gravitational field is reproduced by Eq. (12.121) for two mass points 
is quite similar though the calculation is much more complicated. Again, the contribution 
from the source does not quite reproduce the field satisfying the appropriate gauge condition 
except in the Landau-like gauge (a = 0) though it does satisfy the linearized Einstein 
equations for all a. We verified that the contributions from the source and initial data 
together give the correct field in the causal future of the mass points and that the field at a 
point not in the causal future of the mass points is reproduced by the field from the initial 
data alone. 



26 



Acknowledgments 

We are indebted to Demian Cho, who inspired us to carry out this work and contributed 
to initial discussions. We are also grateful to Emil Mottola for useful correspondence. We 
used Maple 7 in some of our calculations. 

APPENDIX A: A PROOF OF EQ. 022b 

The definition of the delta function (12.61) is equivalent to the property 

d A xy^ f d A x'^g~'B I (x)5 I r 5 A {x,x')A I '{x l ) = [ d A x^B I (x)A I (x) (Al) 



for any smooth compactly-supported functions B I {x) and A I {x). This implies that the delta 
function satisfies 

8 I /5\x',x) = 5 I r 5*{x,x'). (A2) 
Now, for any smooth compactly-supported function Aj(x) we find by Eq. (12.71) 

L T X J J #rfy/=g'G%{x, x')L^, J ' 'Aj,{x') = L J X J 'Aj(x) . (A3) 

Let B ! (x) be an arbitrary smooth compactly-supported function and let fi{x) be the ad- 
vanced solution to the equation L IJ fj = B 1 . Multiplying Eq. (1A3j) by fi(x) and integrating 
over x we have 

Jd A x^9 Jd 4 x'^f I (x)Ll J G« I ,(x,x% i y'Aj,(x') = J d A x^ fi(x)Li J Aj(x) . 

(A4) 

We can integrate by parts in x and x' so that L l J and L l x , J act on fj(x) and Gjj,(x,x'), 
respectively, with the following result: 

Jd A x^B I (x) Li, 3 'G%,(x,x') A P (x') = J d A x^B I (x)A I (x). (A5) 

This equation is equivalent to Eq. (12.71) . 

APPENDIX B: SOME FORMULAS USED IN SEC. [V] 

From Eq. (O]) we find, using Eqs. (10I]) - (O7)) . 

V a (d b cos Hfi) = H\l-2z)g ah} (Bl) 
V a {d v cos Hfi) = H 2 [g ab , + 2(1 - z)n a n b ,} , (B2) 
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which enable us to calculate the derivatives of C^ a , b , and C^, w , i = 1,2,3, defined by 
Eqs. (EMMmU as 



a'b' 



ab 

^'(c^;' b b ' + cK b : 



AH 



a' V 



n c'(a' &M , „ c' fa' 6') 
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x{(2z-l)g c '( a 'n a n b n b '^ 

-g{a c 'n b) n a 'n b ' - 2(1 - z)n a n b n a ' n b V'} 



(B6) 



(B7) 



These are used to find S^ b given by Eqs. fl5.63l) - fl5.67l) . 

For any integrand f(r', cosd') , Eqs. (15.681) - ( T5. 701) can be used to carry out the ^'-integral. 
We define 



Then, 



2ff' 
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(BIO) 
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where we use e = — 1 (e = 1) to calculate the field in the region which is (not) in the causal 
future of the mass point, and 

poo pi 

df / dcos§' f(f', cos 9')9(z- 1) 



if f > A - A' ; 



,f+(A-A') pi 

df \ dcos9' f(f, cos 9'), 
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k yj-rH-(A-A') 

The expression of the form 
evaluated as follows: 
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(B13) 

appearing in these formulas are 



F(f , cos 9')5(1- z) 
/-\ 2AA' 



cos 6'=-l 



A -A' 



x 



5(f' - f - (A - A')) + 5(r -f + (A-A'))-5(r +r-(A-A')) , (B14) 
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and 



F(f f , cos e')5'(l -z) 
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- F ir\ 1)^7^7 ^ + ^ " f + ( A " A '))l 

- F (f , JL f^L 5 (r> + f - (A - AO) 
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